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The recent WMAP and Planck data have confirmed that exotic dark matter together with the 
vacuum energy (cosmological constant) dominate in the flat Universe. Many extensions of the 
standard model provide dark matter candidates, in particular Weakly Interacting Massive Particles 
(WIMPs). Thus the direct dark matter detection is central to particle physics and cosmology. Most 
of the research on this issue has hitherto focused on the detection of the recoiling nucleus. In this 
paper we study transitions to the excited states, possible in some nuclei, which have sufficiently low 
lying excited states. Examples considered previously were the first excited states of and '^^^Xe. 
We examine here ®®Kr, which offers some kinematical advantages and is currently considered as a 
possible target. We find appreciable branching ratios for the inelastic scattering mediated by the 
spin cross sections, with an inelastic event rate of 4.4 x 10“‘^kg“^d~^. So, the extra signature of the 
gamma ray following the de-excitation of these states can, in principle, be exploited experimentally. 
A brief discussion of the experimental feasibility is given. 

PACS numbers: 95.35.+d, 12.60.Jv 11.30Pb 21.60-n 21.60 Cs 21.60 Ev 


I. INTRODUCTION 

The combined MAXIMA-I [l|, BOOMERANG Q, DASI Q and COBE/DMR Cosmic Microwave 
Background (CMB) observations Q imply that the Universe is flat Q and that most of the matter 
in the Universe is Dark [^, i.e. exotic. These results have been confirmed and improved by the 
recent WMAP Q and Planck Q data. Combining the data of these quite precise measurements one 
finds: 


Ufa = 0.0456 ±0.0015, Ocdm = 0.228 ± 0.013, Da = 0.726 ± 0.015 

(the more recent Planck data yield a slightly different combination Dcdm = 0.274 ± 0.020, Da = 
0.686 ± 0.020). It is worth mentioning that both the WMAP and the Plank observations yield 
essentially the same value of but they differ in the value of h, namely h = 0.704 ± 0.013 

(WMAP) and h = 0.673 ± 0.012 (Planckb Since any “invisible” non exotic component cannot 
possibly exceed 40% of the above Dcdm [SI, exotic (non baryonic) matter is required and there 
is room for cold dark matter candidates or WIMPs (Weakly Interacting Massive Particles). Even 
though there exists firm indirect evidence for a halo of dark matter in galaxies from the observed 
rotational curves, see e.g. the review 0: it is essential to directly detect such matter in order to 
unravel the nature of the constituents of dark matter. At present there exist many such candidates: 
the LSP (Lightest Supersymmetric Particle) [IMl, technibaryon [TsL fl^. mirror matter (^. [^. 
Kaluza-Klein models with universal extra dimensions [12, etc. Additional theoretical tools are 
the structure of the nucleus, see e.g. [23 - [27l| . and the nuclear matrix elements [22 - l32 |. 

In most calculations the WIMP is supposed to be the neutralino or LSP (lightest supersymmetric 
particle), which is assumed to be primarily a gaugino, usually a bino. Models which predict a 
substantial fraction of higgsino lead to a relatively large spin induced cross section due to the Z- 
exchange. Such models tend to violate the LSP relic abundance constraint and are not favored. Some 
claims have recently been made, however, to the effect that the WMAP relic abundance constraint 
can be satisfied in the hyperbolic branch of the allowed SUSY parameter space, even though the 
neutralino is then primarily a higgsino . We will not restrict ourselves to supersymmetry and we 
will adopt the optimistic view that the detection rates due to the spin may be large enough to be 
exploited by the experiments, see, e.g., [12, . Such a view is further encouraged by the fact 
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that, unlike the scalar interaction, the axial current allows one to populate excited nuclear states, 
provided that their energies are sufficiently low so that they are accessible by the low energy LSP, a 
prospect proposed long time ago and considered in some detail by Ejiri and collaborators [s^- 
For a Maxwell-Boltzmann (M-B) velocity distribution the average kinetic energy of the WIMP is: 


(T) 


50 keV 


100 GeV 


( 1 ) 


So for sufficiently heavy WIMPs the average energy may be adequate to allow scattering to low 
lying excited states of of certain targets, e.g. of 57.7 keV for the 7/2+ excited state of the 39.6 
keV for the first excited 3/2+ of ^^®Xe, the 35.48 keV for the first excited 3/2+ state of ^^®Te and 
the 9.4 keV for the first excited 7/2+ state of ®^Kr . 

The first requirement is to have a large elementary nucleon cross section. In most models one can 
have both the scalar and spin induced cross section. In the case of the standard recoil experiments 
the coherent cross section is expected to dominate over the spin cross section especially for large 
A, unless the spin induced amplitude is very large at the quark level. This is possible in some 
models like supersymmetry in the co-annihilation region [3^ . where the ratio of the spin to coherent 
nucleon cross section, depending on tan (3 and the WIMP mass, which is in the range 200-500 GeV, 
can be as large as 10^. In a region of the model space the ratio of the elastic spin cross section 
to the coherent can be as large as 10%. More recent calculations in the supersymmetric 5'O(10) 
model [40j |. also in the co-annihilation region, predict large spin to coherent cross section ratio, of 
the order of 2 x 10^ and a WIMP mass of about 850 GeV. Models of exotic WIMPs, like Majorana 
particles with spin 3/2 [4l|, also can lead to large nucleon spin cross sections, which satisfy the relic 
abundance constrain. These large nucleon cross sections may also allow inelastic scattering to the 
excited states with observable rates. In fact we expect that the branching ratio to the excited state 
will be enhanced in the presence of non zero energy threshold, since only the total rate to the ground 
state transition will be affected and reduced, while this will have a negligible effect on the rate due 
to the inelastic scattering. The analysis is simplified if, as expected from particle models as well as 
the structure of the nucleon, the isovector nucleon cross section is dominant. 


II. THE SPIN DEPENDENT WIMP-NUCLEUS SCATTERING 

The spin dependent WIMP-Nucleus cross section is typically expressed in terms of the WIMP- 
nucleon cross section, which contains the elementary particle parameters entering the problem at 
the quark level. For our purposes it is adequate to have at the quark level two parameters, one 
isoscalar 00 ( 9 ) and one isovector ai(g). In going to the nucleon level one must transform these two 
amplitudes by suitable renormalization factors given in terms of the quantities Ag prescribed by 
Fllis HI, namely Au = 0.78 ± 0 . 02 , Ad = -0.48 ± 0.02 and As = -0.15 ± 0 . 02 , i.e. 

oo = ao{q) {Au + Ad + As) = 0.13 ao{q), 

ai = ai{q) {Au — Ad) = 1.26 ai{q). ( 2 ) 

In other words the isovector component is renormalized as is expected for the axial current, while 
the isoscalar component is suppressed, consistent with the FMC effect, i.e. the fact that only a tiny 
fraction of the spin of the nucleon is coming from the spin of the quarks. Thus in general, barring 
very unusual circumstances at the quark level, the isovector component is expected to be dominant. 
It is for this reason that we started our discussion in the isospin basis. From these amplitudes one 
is able to calculate the isoscalar and isovector nucleon cross section. After that, within the context 
of a nuclear model, one can obtain the nuclear matrix element: 

\ME\^ = afSii^u) + aiaoS'oi(u) + alSooiu), (3) 

where Sij are the spin structure functions, which depend on the momentum transferred to the 
nucleus. It is useful to define the structure functions 
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with ill, r^o the isovector and isoscalar static spin matrix elements. Then the functions Fij take 
the value unity at momentum transfer zero. Furthermore these new structure functions are approx¬ 
imately the same [^ . 1^ . Then to a good approximation 


T T ^o^o) 


( 4 ) 


Thus, leaving aside the factor Fn, which will be taken care of explicitly in the expression for the 
event rate (see below), the nuclear static spin cross section becomes 


^spin _ 

~A — 


(fljai^N) + 2sign{aiao)^o^iV^o{N)cri{N) + ngcro(^)) 


( 5 ) 


where a-o{N) and cri(iV) are the elementary nucleon isoscalar and isovector cross sections. Instead 
of these two one could have written an equivalent expression in terms of the proton and neutron 
elementary cross sections [1^ ap and an- 

The above expressions look complicated. In most models, however, e.g. in the case the WIMP is 
the Lightest Supersymmetric Particle (LSP), one encounters Z exchange and Z — qq coupling are 
purely isovector. The same is true in models involving heavy neutrinos as WIMPs, e.g. Kaluza Klein 
theories in models with Universal Extra Dimensions [^. In the latter case it is possible to have both 
an isoscalar and isovector component if the WIMP happens to be a vector boson. Only isovector 
amplitudes are predicted also in the case of exotic spin 3/2 particles [4l|. Since the isoscalar nucleon 
cross section is expected to be negligible, the above expression reduces to 


spin 

~A 


0?ai(iV) 


( 6 ) 


Regarding the spin nucleon cross sections at this point there exist some experimental limits, namely 
for ^^®Xe and ^^^Xe and The first group extracts a limit on the neutron cross section 

of cr„ = 2 X 10“"*^°cm^ = 2 X 10“"*^pb, while from the latter the limit Cp = lxlO“^®cm^ = 1.0xl0“^pb. 
These limits were based on nuclear physics considerations, namely selecting the predominant nuclear 
spin matrix elements. Thus they cannot exclude the presence of the other elementary nucleon cross 
section. On the other hand, if the elementary amplitude is purely isovector, these limits would imply 
CTi = ap+an = 1.02 X 10“^pb. We should mention in passing that the nuclear matrix elements for ^®F 
are expected to be much more reliable [^ . 1^ . Actually this problem will remain open until all the 
three nucleon cross sections (scalar, proton spin and neutron spin) can be determined along the lines 
previously suggested , after sufficient experimental information on at least three suitable targets 
becomes available. Anyway in the present work we will not commit ourselves to any particular model, 
but for orientation purposes we will use the value [4l| of ai{N) fv 1.7 x 10“®®cm“^ = 1.7 x 10“^pb. 


III. THE STRUCTURE OF THE ®®KR NUCLEUS 


A complete calculation of the relevant spin stucture should be done alon g th e lines previously done 
for other targets [l^, [5l|, [l^, ^|, [IJ, and more recently [131,[5^ in the the shell model 

framework. We will not concern ourselves with other simplified models, e.g schemes of deformed 
rotational nuclei [^ . 

A summary of some nuclear ME involved in elastic and inelastic scattering can be found elsewhere 

4il, [13. 

In principle one may have to go a step further in improving the structure functions, since it has 
recently been found that in the evaluation of the ground state (gs) spin structure functions some 
additional input is needed, namely: 


One has to consider the nucleon form factor 


J = 

9A 


|.|43l| entering the isovector axial current: 

((T.q)q 


S = cr- 


q- 


( 7 ) 
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If we choose as a z-axis the direction of momentum we find: 



Thus only the longitudinal component of the transition operator is modified by the inclusion 
of the nucleon form factor. Thus, even at sufficiently high momentum transfers, only 1/3 of 
the differential rate will be affected. 

• One has to consider effects arising from possible exchange currents. 

Such currents lead to effective 2-body contributions obtained in the context of Chiral Effective 
Field Theory (EFT) [Hi. 

We will elaborate further on this point, however, since it has been shown that such effects are 
independent of the target nucleus [H- Thus they can be adequately contained in a retardation 
factor of the isovector amplitude, which takes the value of about 0.8. Thus one can absorb such 
effects in the isovector nucleon spin cross section, which anyway in practice it may have to be 
extracted from experiment, once the spin induced rates have been observed. We expect such a 
treatment will also work in the structure functions involving transitions to excited states. 


IV. SHELL-MODEL INTERPRETATIONS 

To interpret the experimental data, shell-model calculations have been carried out in the 28- 
50 valence shell composed of the lp 3 / 2 , O/ 5 / 2 , lpi/ 2 , O 59/2 orbitals. The calculations have been 
performed with recently available effective interaction by Brown and Lisetskiy [H shell model 
code NuShellX (H- The single-particle energies employed in conjunction with the jj44b interaction 
are -9.6566, -9.2859, -8.2695, and -5.8944 MeV for the IP 3 / 2 , O/ 5 / 2 , lpi/ 2 , 0 ^ 9/2 orbitals, respectively. 
The ground state is correctly reproduce by present calculation. The first excited 7/2+ state is 
predicted at 100 keV, while corresponding experimental value is 9 keV. 

In Table U we have compared experimental B(E2) and i?(Ml) values. The overall agreement is 
reasonable, at least for B(E2) . The calculated il(E2 : 7/2+ —9/2)^) value is 11.62 W.u., while the 
corresponding experimental value is 21.3(17) W.u. The structure of 9/2)^ state is 2 ^( 3972 )”^ (with 
probability ~ 16%). The same configuration is predicted for the first excited state (^ 26%). The 
calculated i?(Ml : 7/2+ —>■ 9/2+) value is 0.0028 W.u., while the corresponding experimental value 
is 0.0095 W.u. The calculated values of the quadrupole and magnetic moments are shown in Table 
im These results are in reasonable agreement with the available experimental data. In particular, 
they reproduce correctly the sign of the quadrupole and magnetic moments. 

Table I: Calculated B{E2) and B{M1) values for ®^Kr isotope with standard effective charges: 

egjj=1.5e, e/g=0.5e and for (the experimental 7 -ray energies corresponding to these 

transitions are also shown). 




B(E2)(W.u.) 


B(Ml)(W.u.) 


Jr ^jj 

Ej (keV) 

Expt. 

jj44b 

Expt. 

jj44b 

7 / 2 / ^ 9 / 2 / 

9.4 

21.3(17) 

11.62 

0.0095 

0.0028 


Table II: Electric quadrupole moments, Qs (in eb) and magnetic moments, ^ (in /tat), (the 
effective charges ep=1.5, e„=0.5 and for gf® = ( 7 ®®®are used in the calculation). 


J" 

Q s,ex.p 

Qs,jj 44 b 

/^exp 

Mjj44b 

■ 97 ^ 

+0.26 (3) 

+0.34 

-0.970669 (3) 

-1.412 

7/2+ 

+0.495 (10) 

+0.41 

-0.943 (2) 

-1.099 
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Figure 1: Elastic spin structure function Flj{Eii) for the isospin modes i,j = 11,01,00 as 
functions of the energy transfer in keV in the case of the target ®^Kr (a). For comparison we 
present the square of the form factor entering the coherent mode (solid line). We also show the 
inelastic spin structure function FI^{Er) (b) . Note that the elastic structure functions, when 
normalized to unity at momentum transfer zero, are essentially independent of isospin. In the case 
of the inelastic scattering the structure functions are not very much depressed as function of the 
energy transfer, compared to those of the elastic scattering. The solid line, dotted line and dashed 
line correspond to ^ 11,^01 and Fqo respectively. 


The thus obtained static spin MEs are 

rio = 1.037, fli = —1.018 for elastic transitions, 

flo = —4.880 X 10“^, fli = 4.421 x 10“^ for inelastic transitions. 

Looking at these matrix elements we notice the following: 

• The nuclear wave functions do not favor the isoscalar contribution to overcome the expected 
suppression of the corresponding isoscalar amplitude in going from the quark to the nucleon 
level discussed above. 

• The inelastic transition is not predominantly a Gamow-Teller-like transition (that is, driven 
by the a operator), but but it is dominated by higher multipoles. 

The disadvantage of small value of the spin ME in the inelastic transition may, however, be 
partly overcome by the fact that, unlike the EM decay of the excited state, a relatively high 
momentum can be transferred to the nucleus by the WIMP and the excited state structure 
function does not fall fast with the transferred momentum. 

The structure functions relevant to our calculation are exhibited in Fig. [T] 


V. THE FORMALISM FOR THE WIMP-NUCLEUS DIFFERENTIAL EVENT RATE 


The expression for the elastic differential event rate is well known, see e.g. [^,[431 . For the 


reader’s convenience and to establish our notation we will briefly present the essential ingredients 
here. We will begin with the more familiar time averaged coherent rate, which can be cast in the 
form: 


di?o 

dEfi 


A 


Px 

Arup 



V< > 


1 







F'^{u) 4'o(aVM) 


( 9 ) 
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with with ilJ-p) the WIMP-nucleus (nucleon) reduced mass and A is the nuclear mass number. 

is the WIMP mass, p(x) is the WIMP density in our vicinity, assumed to be 0.3 GeV cm“^, 
and rrit the mass of the target, u is the recoil energy En in dimensionless units introduced here 
for convenience, u = = AmpEn, with A the nuclear mass number of the target and b the 

nuclear harmonic oscillator size parameter characterizing the nuclear wave function. It simplifies 
the expressions for the nuclear form factor and structure functions. In fact: 


Er 

QoiA) 


QoiA) = [nipAb^]-^ = 40y4-4/3 MeV 


( 10 ) 


The factor in the above expression is vo /since in Eq. ([9]) the WIMP flux is given 

in units of (u^) appears. In the above expressions a = {y/2firbvo)~^, vq the velocity of the sun 
around the center of the galaxy and F(u) is the nuclear form factor. Note that the parameter a 
depends both on the WIMP mass, the target and the velocity distribution. 

For the axial current (spin induced) contribution one finds for the elastic WIMP-nucleus scattering: 


dRp _ mt 

(IEr ^ Arrip 



V< > 


1 


a 


spin 

A 



Spin 




Eii(u) ’i>o(a^/u) 


( 11 ) 

where En is the isovector spin response function, i.e. Egg for the elastic case and Fin for the inelastic 
one. Their behavior is exhibited in Fig. [1] 

We notice that the only difference in the shape of the spectrum between the coherent and the spin 
comes from nuclear physics. Since the square of the form factor and the spin structure functions 
have approximately similar shapes, we expect the corresponding differential shapes to be the same, 
and only the scale to be different. 

The function 4'o(a:) is defined by 'I'o(a:) = g(vmin,VE{oi))/vp, where vp is the velocity of the sun 
around the center of the galaxy. g{vmin,'UE{oi)) depends on the velocity distribution in the local 
frame through the minimum WIMP velocity for a given energy transfer, i.e. 


'^min 


ArupER 


2 p 


2 

r 


( 12 ) 


The above way of writing the differential event rates we have explicitly separated the three im¬ 
portant factors: 


• the kinematics, 

• the nuclear cross section or crj‘" 

• the combined effect of the folding of the velocity distribution and the form factor or the nuclear 
structure function. 


For the Maxwell-Boltzmann (M-B ) distribution in the local frame it is defined as follows: 

1 f 


gig^min^^E^ — 


'^max — '^esc'i 


{y/^vp)" 




(13) 


where is the velocity of the Earth, including the velocity of the sun around the galaxy. We have 
neglected the velocity of the Earth around the sun, since we ignore the time dependence (modulation) 
of the rates. The above upper cut-off value in the M-B is usually put in by hand. Such a cut-off 
comes in naturally, however, in the case of velocity distributions obtained from the halo WIMP mass 


density in the Eddington approach [63|, which, in certain models, resemble a M-B distribution [64 1. 


Integrating the above differential rates we obtain the total rate including the time averaged rate 
for each mode given by: 


D _ 

-*^coh 


Px rnt 


Amp 


\/< > A^ E 


Ieoh 


r(yesc/t^) 

lEth/QoiA) 


dt 

du 


du, 


(14) 
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R. 


Px 

m-)^ Arup 



2 

\/< 


r (Vesc /fi) 
JEth/Qo{A) 


dt 

du 


du 


spin 


(15) 


for each mode (spin and coherent). Etui A) is the energy threshold imposed by the detector. 

These expressions contain the following parts: i) The gross properties and kinematics ii) The 
parameter t, which contains the effect of the velocity distribution and the nuclear form factors 
iii) The WIMP-nuclear cross sections A^cfn or crj‘”. The latter, contains the nuclear static spin 
ME. From the latter the elementary nucleon cross sections can be obtained, if one mode becomes 
dominant as already mentioned above. Anyway using the values for nucleon cross sections, in 
Eq. (full and crj”’ in Eq. (ITS)) , we can obtain the total rates. 

Conversely, if only one mode is dominant, one can extract from the data the relevant nucleon cross 
section (coherent, spin isoscalar or spin isovector) or obtain exclusion plots on them. 


VI. TRANSITIONS TO EXCITED STATES 

Transitions to excited states are normally energetically suppressed, except in some odd nuclei that 
have low lying excited states. Then spin mediated transitions to excited states are possible. The 
most favorable are expected to be those that, based on the total angular momentum and parity of 
the states involved, appear to be Gamow-Teller like transitions. A good possibility is the 7/2+ state 
of ®^Kr, which is at 9.5 keV above the 9/2+ ground state. 


A. Isotope considerations 

Transitions to the first excited states by spin interacting WIMPs (SD WIMPs) can occur in the 
case of some odd A targets, if the relevant excitation energy is E < 100 keV. This is due to the 
high velocity tail of the M-B distribution, so that a reasonable amount of the WIMP energy may be 
transferred to the recoiling nucleus. 

Possible odd mass nuclei involved in DM detectors used for WIMP searches are the 57.6 keV 
state in ^^^I and the 39.6 keV state in ^^®Xe. The spin excitations of these states are not favored, 
because the dominant components of the relevant wave functions are characterized by A£ ^ 0 , i.e. 
£ forbidden transitions. Nevertheless the spin transitions are possible, due to the small components 
as seen from the Ml 7 transition rates. Other possibilities are, of course, of interest, e.g. ®^Kr is a 
good experimental candidate. One can easily enrich it isotopically and it can make a good liquid or 
gas detector. Thus it is quite realistic to study the inelastic excitations in this nucleus to search for 
SD WIMPs. 

In fact the experimental observation of the inelastic excitation has several advantages. The ex¬ 
perimental feasibility is discussed in section VIII. 

Table III: Inelastic spin excitations of ®^Kr. A : natural abundance ratio, E : excitation energy, 

ground state spin parity, J/ : excited state spin parity, and Ti/ 2 '. half life. For comparison we 

present the same quantities for ^^^Te. 


Isotope A(%) E (keV) 

Ji Jf 

ri/2(ns) 

*'=‘Kr 11.5 

125 Te 7.07 

9.5 

35.5 

9/2+ 7/2+ 
1/2+ 3/2+ 

154 

1.48 








B. Kinematics 


The evaluation of the differential rate for the inelastic transition proceeds in a fashion similar to 
that of the elastic case discussed above, except: 

1. The transition spin matrix element must be used. 

2. The transition spin response function must be used. For Gammow-Teller like transitions, it 
does not vanish for zero energy transfer. So it can be normalized to unity, if the static spin 
value is taken out of the ME. 


3. The kinematics is modified. The energy-momentum conservation reads: 


+ v^q — Ex=0, Fix = excitation energy - Er + v^\J 2mAEfi — = 0 (16) 


Clearly ^ > 0 as before. Then ^ < 1 and the reality of Eji impose the condition: 


V > 


Ex + ^Er 


(17) 


We thus find that for a given energy transfer En the minimum allowed WIMP velocity is given 
by: 


Ex + ^Er 

_ _Mr_ 

^2mAER 


(18) 


while the maximum and minimum energy transfers are limited by the escape velocity in the 
WIMP velocity distribution. We find that: 


(.^i?)mm ^ Er ^ {^Er) max 


(19) 


with 


Mr ( 2 Ex 


niA 


fir 


Ex 

/ir 


= — I - — -\ vte- 2rLc— , (Er)^^^ = el. ^ 


Mr ( 2 Ex 


Ex 


niA 


( 20 ) 


In the case of elastic scattering we recover the familiar formulas: 

{ER^rnin — O5 iER)max — 2 '^esc 

ruA 

From the above expressions it is clear that for a given nucleus and excitation energy only 
WIMPs with a mass above a certain limit are capable of causing the inelastic transition, i.e. 


. I 2 472^4 , 1 

m^>mA{ 2^esc^ - 1 I 


( 21 ) 


(m^) . = (4.6, 19, 34, 21) GeV for “'^Kr, respectively 


We find it simpler to deal with the phase space in dimensionless units. Noticing that u = 
{l/2)q%‘^ and 




( 22 ) 
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we find: 


J q^d^dqS + v^q - 

i.e. we recover the same expression as in the case of ground state transitions. 
The above constraints now read: 


y> a 


U + Uo 



(23) 


(24) 


^0 — 


\/2lIrVob’ 


Er 

Q^y 


(25) 


It should be stressed that for transitions to excited states the energy of recoiling nucleus must 
be above a minimum energy, which depends on the escape velocity, the excitation energy and 
the mass of the nucleus as well as the WIMP mass. This limits the inelastic scattering only 
for recoiling energies above the values {Er)^^^. The minimum and maximum energy that can 
be transferred is: 


^min 


1 

4 




(26) 


The maximum energy transfers Umax depend on the escape velocity Vexc = 620km/s. Here we have 
denoted Vexc = VescVQ, where vg is the characteristic velocity of the M-B distribution taken to be 
220 km/s, i.e. yesc ~ 2.84. The values {ER)min and (Er)^^^ relevant for the inelastic 

scattering of ®^Kr and ^^^Te are shown in Table IlYl The minimum required and the maximum 
allowed energy transfer depend on the WIMP escape velocity. The dependence of the maximum is 
not crucial, since the contribution to the total rate becomes negligible at high energy transfers. The 
behavior of the minimum in the range of the existing limits of the escape velocities is exhibited 
in Fig. [2] 


Table IV: The kinematical parameters entering the inelastic scattering to the first excited state of 
®^Kr. For comparison we present the same quantities for ^^^Te. 


target 

parameter 

m;,(GeV) 



20 

50 

100 

200 

500 

1000 



0.671 

0.645 

0.638 

0.633 

0.632 

0.631 


Umax 

0.440 

1.71 

3.67 

6.17 

9.12 

10.5 


(^^«)min(keV) 

0.078 

0.049 

0.049 

0.049 

0.048 

0.048 


(-E«)ma,x(keV) 

51.6 

199 

425 

716 

1059 

1228 

^'"Te 

^min 

0.014 

0.012 

0.012 

0.011 

0.011 

0 .0.11 


Umax 

0.45 

2.35 

8.95 

11.69 

20.00 

24.86 


(■E«)min(keV) 

0.94 

0.78 

0.74 

0.72 

0.71 

0.71 


(^lfl)Tna,x(keV) 

30.2 

157 

400 

786 

1347 

1670 


To avoid uncertainties arising from the relevant particle model, we will present the rate to the 
excited relative to that to the ground state (branching ratio). The differential event rate for inelastic 
scattering takes a form similar to the one given by Eq. [TT] except that 


u + uo\ 

^7^)' 


Hi ^ Fii{u) -)■ Fii(«)”'*'““% 4'o(aVw) ^ ^'o 
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Figure 2: The minimum required energy transfer as a function of the escape velocity for ®^Kr (a). 
For comparison we present the same quantity for ^^®Te (b). 


VII. SOME RESULTS 

For purposes of illustration we will employ the nucleon cross section of 1.7 x 10“^pb obtained in a 
recent work , without committing ourselves to this or any other particular model. Another input 
is the WIMP density in our vicinity, which will be taken to be 0.3 GeV cm“^. Finally the velocity 
distribution with respect to the galactic center will be assumed to be a M-B with a characteristic 
velocity vq = 220km/s and an upper cut off (escape velocity) of 2.84r;o. 
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Er —^keV 

Figure 3: Energy spectrum for WIMP ®^Kr elastic scattering (a) and that for the 9.4 keV excited 
state inelastic scattering (b). The thick solid, solid, dotted, dashed-dotted and dashed curves 
correspond to WIMP masses 50, 100, 200, 500 and 1000 GeV. 



(a) 



Er —^keV 


(b) 


Figure 4: The ratio of the differential scattering rates, j of 10 as a 

function of the recoil energy Er in keV in the case of the target ®^Kr. The thick solid, solid, 
dotted, dashed-dotted and dashed curves correspond to WIMP masses 50, 100, 200, 500 and 1000 
GeV. In panel (a) we show the results for j/esc = 2.5 and in (b) for j/esc = 2.8. The range of the 
allowed energy transfer, depends, of course, on the WIMP mass. The dependence on the escape 
velocity in the range of the accepted values is mild. 


A. The differential event rates 

The differential event rates, perhaps the most interesting from an experimental point of view, 
depend on the WIMP mass, but we can only present them for some select masses. Our results for 
the elastic differential rates for typical WIMP masses are exhibited in Fig. [5^. For comparison we 
present the differential event rates for transition to the excited state in Fig. [3 ]d. 

We will evaluate the branching ratio of the inelastic differential cross section to that of the ground 
state ignoring the coherent mode. We will restrict ourselves in the isovector transition. This is 
reasonable, since as we have already mentioned the isoscalar is absent in most particle models and 
even if it appears at the quark level it is expected to be suppressed due to considerations related 
to the spin of the nucleon. In such a case the ratio becomes independent of the elementary nucleon 
cross section. The obtained results are exhibited in Fig. S) 
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Figure 5: The ratio of the total scattering rates, i?(excited)/i?(gs), as a function of the WIMP 
mass in GEV for y = 2.84 in the case ®^Kr. From bottom up the threshold values are 0,1, 2, 3,4, 7 
and 10 keV. Only the spin mode has been taken into account. No quenching factor has been used. 


fcuO 


O 



^WIMP — ^G6V 


Figure 6: We show for the target ®^Kr the time average total rate as a function of the WIMP 
mass in GeV obtained in the case of elastic scattering (a) and inelastic scattering (b), assuming a 
zero energy threshold. Note that the maximum in the case of the inelastic scattering occurs at a 
higher mass compared to the location of the maximum for the ground state transition. 


B. Total rates 

From expressions (fHl) and csi), we can obtain the total rates. It is instructive to start with a 
consideration of the branching ratio of the total rates, which for the reasons discussed above is going 
to be independent of the elementary nucleon cross section. Since the elastic scattering event rate 
is reduced by the threshold effects, but the inelastic scattering is not affected by such effects, we 
expect the branching ratio to be increasing as the threshold energy is increasing. The situation is 
exhibited in Fig. [S] 

The total rates obtained assuming zero energy threshold are exhibited in Fig. [6] as functions of 
the WIMP mass. 
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From this figure we see that the maximum total rate for the inelastic transition is around a WIMP 
mass of about 50 GeV, while the location of the maximum for ground state transition is around 20 
GeV. On the other hand the excited state cannot be populated in the case of WIMPs lighter than 
4.6 GeV, but the ground state can, albeit with small rate there. We also see that for a WIMP mass 
of 20 GeV the obtained rate is 0.16 counts/ (kg-Y')= 4.4 x 10“"^/(kg-d). This is just above the rate of 
2.8 X 10“"^ counts/(kg-d) reported previously . This is a little surprising because a larger spin 
\ME\^ with the value of 0.8 x (0.08/1.91^) = 0.0175 was estimated in the earlier work. The larger 
rate obtained in the present work could be attributed to our use of an energy dependent structure 
function, or, most likely, to the supersymmetric model they used, which would have yielded a nucleon 
cross section smaller than the one used here. 


VIII. EXPERIMENTAL ASPECTS OF INELASTIC NUCLEAR SCATTERING RATES 

In this section, we discuss experimental aspects of SD WIMP studies by measuring inelastic 
nuclear scattering. So far SD and SID WIMPs have been studied experimentally by measuring 
nuclear recoils of elastic scatterings. 

SD WIMPs may show fairly appreciable cross sections of inelastic spin excitations, as shown in 
previous sections. Experimentally, inelastic nuclear excitations provide unique opportunities for 
studying SD WIMPs. Inelastic excitations are studied in two ways: singles measurement of both the 
nuclear recoil energy and the decaying y-ray energy E^ in one detector, and coincidence measure¬ 
ment of the nuclear recoil and the y-ray in two separate detectors. This is done in a fashion analogous 
to that of the earlier analysis [13 ■ So far SD and SID WIMPs have been studied experimentally by 
measuring nuclear recoils of elastic scatterings. Experimentally, inelastic nuclear excitations provide 
unique opportunities for studying SD rates for WIMP-nuclear interactions. Inelastic excitations are 
studied by two ways. A: singles measurement of both the nuclear recoil energy E^ and the decaying 
y-ray energy E^ in one detector, and B: coincidence measurement of the nuclear recoil and the y-ray 
in two separate detectors. The merits of each of them are as follows. 

The large energy signal is obtained by summing the nuclear recoil signal and the y ray signal. It 
is given as 


E{ex) = E^ + QiERiex))En{,ex), (27) 

where E/j(ex) is the nuclear recoil energy, E^ is the excitation energy and Q{ER{ex)) is the quenching 
factor for the recoil energy signal. In most scintillation and ionization detectors, the quenching 
factor is as small as Q{ER{ex)) Ri 0.1 — 0.05. It must be determined for each target and detector 
experimentally. Eor ®^Kr, the primary target considered in this work, the quenching factor is not 
important for a gaseous detector, but for a liquid or solid it is taken to be 0.08 — 0.1. An overall 
picture can be obtained by a phenomenological approach based on the Lindhard theory , 

which is exhibited in Pig. 0 In the case of ®^Kr, due to the low excitation energy, our conclusions 
about the inelastic transition are not affected much by the quenching factor. In the present work 
we will use a quenching factor of 0 0.08. 

Therefore it appears that the energy deposited is mainly the excitation energy. This is much 
larger than just the recoil energy signal of E{gi) = Q{ER{gT)), which is quenched, depending on the 
detector. 

The sharp rise of the energy spectrum at the energy of E^ + Q{Emin)Emim where Emin is the 
minimum energy transfer to the recoil nucleus. This makes it possible to identify the WIMP nuclear 
interaction. On the other hand, the recoil energy spectrum i?/j(gr) is continuum like back ground 
at the low energy region, and thus is hard to be identified. 

E{ex) is well above the detector threshold E{th), while the main part of E{gi) is cutoff by E{th). 
Accordingly, the event rate i?(ex) is about the same order of magnitude as R{gr) for SD WIMPs 
although the inelastic cross section is much smaller than the elastic one. 

The typical energy spectra to be measured experimentally for the elastic and inelastic transitions 
of ®^Kr are shown in Pig. [51 Here we assumed detectors with the quenching factor of (5=0.05 and 
the energy threshold of E{th) = 1.6 keV. The yield on the y axis is the one per unit energy of the 
electron equivalent energy, i.e. QEr and the energy on the x axis is the electron equivalent energy. 
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Figure 7: The quenching factor as a function of the recoil energy in keV obtained in the Lindhard 
theory is presented at low energy (a) and at somewhat larger energies (b). 



(a) (b) 



(c) 


Eft ^keVee 


(d) 


Figure 8: Energy spectrum for WIMP ®^Kr elastic scattering (panels (a) and (b)) in two energy 
regions. Note that the maximum allowed energy transfer depends on the WIMP mass. For a heavy 
WIMP it can reach about 60 keVee, but it falls real fast with the energy transfer. The energy 
spectrum for the 9.4 keV excited state inelastic scattering is also shown (panel (c)). The quenching 
factor employed was (3=0.08 and the energy threshold E{th) = 1.6 keV. The x and y scales are the 
electron-equivalent energy and the rate per unit electron equivalent energy. 
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Er ^keVee 


Figure 9: We show the energy spectrum for the elastic transition (a) and the inelastic to the first 

excited state (b) for the target ^^®Xe 



Er ^keVee 

Figure 10: The same as in Fig. IHlfor the target 


We note that the yield is enlarged by a factor 1/Q=20, while the energy is reduced by a factor 
(5=0.05. The low energy part of the elastic scatting is cut off by the threshold energy of 1.6 keV 
electron equivalent energy, i.e. 32 keV recoil energy. 

For comparison we plot the same data for Xe ^^®Xe and in Figs 1^ and ITUl respectively , which 
have appeared elsewhere 

Comparing the two targets we see that the smallness of the spin ME in the case of ®^Kr is 
compensated by the favorable energy dependence of its inelastic spin structure function and the low 
excitation energy. 


IX. EXPERIMENTAL FEASIBILITY 

The main purpose of the research described in this article is to explore the feasibility of an ex¬ 
periment to search for Cold Dark Matter via an inelastic excitation of a nuclear target. The isotope 
®^Kr was chosen for two reasons: first, it has a low energy Ml transition from the 7/2+ first excited 
state at 9.4-keV to the 9/2+ ground state; second, in the liquid state it is a fairly good scintillation 
detector material. The disadvantage of this choice is that the natural abundance of ®^Kr is only 
11.5% which will require isotopic enrichment. In addition any content of the ubiquitous radioactive 
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isotope ®®Kr will render a sensitive search ineffective. Even the isotopic enrichment of Kr gas from 
the usual sources will never be free enough from ®®Kr to allow a sensitive experiment. Fortunately, 
a recent discovery by other researchers can very probably be used to ameliorate this difficulty. First, 
we discuss the scintillation properties of liquid Kr, and then possible new sources of low-radioactive 
Kr gas. The scintillation and other physical properties of liquid Kr has been well studied and pub¬ 
lished in the literature [i^- [t^. It has been shown by independent groups [13 [nl that the addition 
of a few percent by mass of Xe enhances the fast component of scintillation of LKr by a factor of 
approximately 10. Considering the long lifetime of the 7/2+ to 9/2+ ground-state decay in ®®Kr, 
the pulse sha pe d istortion might be utilized to provide partial background rejection as discussed in 
another case [^. Nevertheless, this advantage would be minimal in the case that the background 
from ®®Kr was not negligible. In the usual case, sources of krypton and argon for industrial use 
are obtained by distillation of atmospheric air, which contains radioactive ®®Kr (10.75y) and ®®Ar 
(269y). These contaminations are unacceptable for use in ultra-low background experiments. It is 
a well-known fact in the geology and geochemistry community that through petroleum exploration, 
very large deep underground accumulations of CO 2 have been discovered in the greater Colorado 
Plateau and in the Southern Rocky Mountain Region 73|. Recently the Princeton-lead Darkside 


Collaboration has been successful in working with industry to extract ultra-low radioactive back- 


m 


ground argon from large volumes of CO 2 using a Vacuum Pressure Swing Absorption plant 
In the case of Kr, this would be the input gas for the isotopic enrichment. These facts imply the 
possibility of a sensitive search for the Cold Dark Matter excitation of ®®Kr to the first excited 
state and the detection of the deexcitation 9.4-keV gamma ray. However, the extraction of the Kr 
from the deep underground CO 2 , as well as the requirement of isotopic enrichment would be far 
more costly than the extraction of ^°Ar, which has an isotopic abundance of 100%. Accordingly, a 
research and development study would have to be undertaken to determine the technical and cost 
feasibility of this experimental approach. However, of the 30 isotopes discussed by Ellis, Flores and 
Lewin ®®Kr is one of the few for which there are clear paths to building a large detector. 


X. CONCLUDING REMARKS 

SI and SD WIMPs have extensively been studied, so far, by measuring elastic nuclear recoils. 
The elastic scattering of SI WIMPs is coherent scatting, thus the cross section is enhanced by the 
factor with A being the nuclear mass number. On the other hand the elastic cross section of 
SD WIMPs is, is in general, smaller by 2-3 orders of magnitude than that for SI WIMPs because 
the spin induced rates do not depend on A^, i.e. they do not exhibit coherence. It may, however, 
compete with the coherent scattering in models in which the spin induced nucleon cross section is 
much larger than the one due to a scalar interaction. We have seen that there exist viable such 
particle models. In such cases the inelastic WIMP-nucleus scattering becomes important. 

Indeed the inelastic scattering via spin interaction provides a new opportunity for studying SD 
WIMPs. Experimentally, observation of both the nuclear recoils and the 7 ray following the excited 
state does lead to the large energy signal of the unquenched and the sharp rise of the energy 
spectrum at around E^. Even though the SD inelastic cross section is smaller than the SD elastic one, 
the inelastic event rate is comparable with the elastic one, since the inelastic signal is well beyond 
the detector threshold energy, while the elastic signal is mostly cut off by the detector threshold. 

In the present paper we discussed mainly the inelastic excitations of ®®Kr. Another possible 
isotope, in addition to the ^^^I and ^^®Xe discussed previously [ 63 , is the ^®Ge, in high energy 
resolution Ge detectors, and ^^®Te currently under study. In short, the present paper, in conjunction 
with the earlier calculations [^ . indicates that the inelastic scatting opens a new powerful way to 
search for SD WIMPs. 
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